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Discovery of an X-ray Nebula associated with PSR 
J2124-3358 
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Abstract. We report the discovery of an X-ray neb- 
ula associated with the nearby millisecond pulsar PSR 
J2124— 3358. This is the first time that extended emis- 
sion from a solitary millisecond pulsar is detected. The 
emission extends from the pulsar to the northwest by 
^0.5 arcmin. The spectrum of the nebular emission can 
be modeled by a power law spectrum with photon index of 
2.2 ±0.4. This is inline with the emission being originated 
from accelerated particles in the post shock flow. 



1. Introduction 

Rotation-powered pulsars are well-known to be the ob- 
jects tapping their rotational energy into pulsed emission. 
However, it is now generally believed that a much larger 
fraction of the rotational energy leaves the pulsars' mag- 
nctosphere in the form of magnetized wind. When the 
relativistic wind particles interact with the interstellar 
medium, synchrotron emission, which is characterized by 
a power-law spectrum (Chevalier 2000) is radiated. 

If a pulsar moves across the surrounding medium with 
a velocity exceeding the speed of sound for the medium, 
a bow shock can be formed. The structure of a bow shock 
can consist of two parts (Gaensler 2005). The forward 
shock, resulting from collisional excitation, is expected to 
produce Ha emission. Also, the termination shock is pro- 
duced by the pressure difference between regions ahead 
of and behind the pulsar's velocity vector. The relativis- 
tic wind particles will be accelerated in the termination 
shock and hence generate synchrotron emission from ra- 
dio to X-ray. The morphology of the X-ray termination 
shock is generally cometary. 

Recently, we have conveyed searches for diffuse X-ray 
emission around a group of millisecond pulsars. The pul- 
sars are PSR J2124-3358, J0437-4715, J0030+0451 and 
J1024— 0719 which have comparable spin parameters (Hui 
& Becker 2006). The period and period derivative of this 
group have ranges of 4.87-5.76 ms and (1.0-1.87)xl0~ 20 
s s _1 respectively In this work, we have discovered an 
elongated structure associated with PSR J2124— 3358 
(Figure 1). 



PSR J2124-3358 was discovered by Bailes et al. (1997) 
during the Parkes 436 MHz survey of the southern sky. 
The pulsar has a rotation period of P = 4.93 ms and a 
proper motion corrected period derivative of P = 1.33 x 
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s s . These spin parameters imply a characteristic 



age of 5.86 x 10 yrs and a dipole surface magnetic field 
of 2.60 x 10 8 G (Manchester et al. 2005). The radio dis- 
persion measure gives a distance of about 250 kpc. It is 
found to have a space velocity of ~ 58 km/s (Manchester 
et al. 2005). Gaensler, Jones & Stappers (2002) discov- 
ered an Ha-emitting bow shock nebula around PSR J2124- 
3358. This bow shock is very broad and highly asymmetric 
about the direction of the pulsar's proper motion. 



2. Data Analysis and Results 

PSR J2124-3358 has been observed by XMM-Newton in 
2002 April 14—15 with an effective exposure time of ~ 40 
ks. A vignetting corrected image of the field of PSR J2124- 
3358 as seen by the XMM-Newton's MOS1/2 CCDs is 
shown in Figure la. The binning factor in this image is 
6 arcsec. Adaptive smoothing with a Gaussian kernel of 
a < 4 pixels has been applied to the image. X-ray contours 
are calculated and overlaid on the image. The contour lines 
arc at the levels of (4.2, 5.2, 7.6, 13, 28, 63) x 10~ 6 cts s" 1 
arcsec -2 . It can be seen that the X-ray source which is co- 
incident with the pulsar position has an asymmetric source 
structure of ~ 0.5 arcmin extent, with its orientation to 
the northwest. Systematic effects which could cause an ad- 
equate distortion of the instrument's point spread function 
(PSF) are not known for XMM-Newton. We are therefore 
prompted to interpret this elongated structure in terms of 
a pulsar X-ray trail. The signal-to-noise of this elongated 
feature in the XMM-Newton data is ~ 4 in the energy 
range 0.25 — 5 keV. 

We have investigated a possible contribution to the 
diffuse X-ray emission by nearby stars. To do so, we in- 
vestigated the Digitized Sky Survey data (DSS) for the 
sky region around PSR J2124— 3358. There are four field 
stars, which are labeled as A, B, C and D in Figure la, in 
the 1.5 arcmin neighborhood of PSR J2124— 3358. None 
of them is found to match the position of the diffuse elon- 
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Fig. 1. a) XMM-Newton MOS1/2 image of PSR 
J2124— 3358 with overlaid contours. The pulsar proper 
motion is indicated by an arrow. The position of bright 
stars located in the 1.5 arcmin neighborhood is indicated, 
b) PSR J2124-3358 as seen by the ACIS detector aboard 
Chandra. 

gated X-ray structure. It can be seen in Figure la that the 
positions of stars B and C coincide with two faint X-ray 
sources which are disconnected with the trail emission of 
PSR J2124-3358, though. 

The detection is which supported by the Chandra 
ACIS-S3 observation took place on 2004 December 19- 
20 for an exposure time of ~ 30 ks. An image made from 
this data with 0.5 arcsec binning and adaptive smooth- 
ing applied (using a Gaussian kernel with a < 1.5 pixel) 
is shown in Figure lb. Arc-like diffuse emission which is 
within the pulsar's H Q nebula is clearly detected. Since 
the PSF width of XMM-Newton is about 10 times that 
of Chandra, it blurred most of the detailed structure seen 
in the Chandra data. However, it should be noted that 
the overall direction of the feature in the Chandra image 
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Fig. 2. Energy spectrum of the X-ray trail of PSR 
J2124-3358 as observed with MOS1/2 and ACIS-S3 de- 
tectors and simultaneously fitted to an absorbed power- 
law model (upper panel) and contribution to the \ 2 R t 
statistic (lower panel). 

is consistent with the orientation of the trail detected by 
XMM-Newton. The signal-to-noisc of this feature in the 
Chandra data is ~ 5 in the energy range of 0.3 — 8 keV. 

For the spectral analysis, we extracted the energy spec- 
trum within a 30 x 35 arcsec box from the XMM-Newton 
MOS1/2 data. Using a box rather than a circular selec- 
tion region allows to avoid the emission from the pulsar 
and excludes any potential contamination from the field 
stars B and C. However, we estimate that the wing of the 
XMM-Newton PSF centered at the pulsar position still 
contributes ~ 20% to the total counts inside the box. The 
background spectrum was extracted from a source free re- 
gion near to PSR J2124-3358. In total, 92 and 67 source 
counts are available from the trail in the MOS1 and MOS2 
detectors, respectively. Response files were computed by 
using the XMMSAS tasks RMFGEN and ARFGEN. The 
spectra were dynamically binned so as to have at least 10 
counts per bin. 

In the Chandra data we selected the energy spectrum 
of the diffuse emission from a box of size 10 x 20 arc- 
sec. Owing to the narrow PSF of Chandra the contamina- 
tion of pulsar emission in this box is negligible. The back- 
ground spectrum was extracted from a low count region 
near to the diffuse feature. In total 46 source counts are 
contributed from the Chandra data. Response files were 
computed by using tools MKRMF and MKARF of CI AO. 
The spectrum was binned to have at least 8 counts per 
bin. The degradation of quantum efficiency of ACIS was 
corrected by applying the XSPEC model ACISABS. 

We hypothesize that the diffused emission should be 
originated from the interaction of pulsar wind and the 
ISM. Synchrotron radiation from the ultra-relativistic 
electrons is generally believed to be the emission mech- 
anism of the pulsar wind nebula, which then is character- 
ized by a power-law spectrum. With a view to test this 
hypothesis, we fitted an absorbed power-law model to the 
nebular spectrum with XSPEC 11.3.1 in the 0.25 - 10 
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keV energy range. With a column density of 5 x 10 20 
cm~ 2 as obtained from spectral fits to the pulsar emis- 
sion, we found that the model describes the observed spec- 
trum fairly well (\l = 0.79 for 26 D.O.F.). The photon 
index is a = 2.2 ± 0.4 and the normalization at 1 keV 
is (2.94 ± 0.48) x 10~ 6 photons keV" 1 cm~ 2 s" 1 (1 - a 
error for 1 parameter in interest). In view of the low 
photon statistic we tested for a possible dependence of 
the fitted model parameters against the background spec- 
trum. All deviations found in repeating the fits with dif- 
ferent background spectra were within the 1 — a confi- 
dence interval quoted above. The unabsorbed fluxes and 
luminosities deduced for the best fit model parameters 
and the energy ranges 0.1—2.4 keV and 0.5—10 keV are 
fx = 1.8 x 10~ 14 ergs s" 1 cm~ 2 , L x = 1.3 x 10 29 ergs 
s- 1 and fx = 1.2 x 10~ 14 ergs s" 1 cm" 2 , L x = 8.9 x 10 28 
ergs s _1 , respectively. The best fitting spectral model is 
displayed in Figure 2. 

3. Discussions 

Adopting the dispersion measure based distance of ~ 250 
pc, the X-ray trail has a length of I ~ 1.1 x 10 17 cm. For the 
pulsar's proper motion velocity of 58 km s _1 (Manchester 
et al. 2005), the timescale, ifl ow , for the passage of the 
pulsar over the length of its X-ray trail is estimated to be ~ 
600 yrs. According to the discussion in Becker et al. (2005) 
on the trail emission of PSR B1929+10 we estimate the 
magnetic field in the shocked region by assuming tn ow to 
be comparable to the electron lifetime of the synchrotron 
emission. This yields ~ 30^G for the inferred magnetic 
field strength in the emitting region. The magnetic field 
strength in the ISM is estimated to be ~ 2 — 6/iG (cf. Beck 
et al. 2003 and references therein). Taking into account 
that the magnetic field in the termination shock might be 
compressed (e.g. Kennel & Coroniti 1984), our estimation 
is approximately consistent if the compression factor is 
- 7. 

Following Becker et al. (2005), we applied a simple one 
zone model (Chevalier 2000) to estimate the spectral be- 
havior and the X-ray luminosity of the nebular emission. 
The X-ray luminosity and spectral index depend on the 
inequality between the characteristic observed frequency 
v x and the electron synchrotron cooling frequency v c 
which is estimated to be 1.6 x 10 17 Hz. Since in gen- 
eral Vx > v c , we concluded that the emission is in a 
fast cooling regime. Electrons with the energy distribu- 
tion, JV(7) oc 7 _p , are able to radiate their energy in the 
trail with photon index a = (p + 2)/2. The index p due 
to shock acceleration typically lies between 2 and 3 (cf. 
Cheng, Taam, & Wang 2004 and references therein). Tak- 
ing p = 2.35 yields a th ~ 2.2 which is in accordance with 
the result from the observed value a obs = 2.2 ± 0.4. As- 
suming the energy equipartition between the electron and 
proton (Cheng, Taam, & Wang 2004), we take the frac- 
tional energy density of electron e to be ~ 0.5 and the 



fractional energy density of the magnetic field es to be 
~ 0.01. Assuming a number density of ISM to be 1 cm~ 3 , 
the distance of the shock from the pulsar is estimated to 
be ~ 3.6 x 10 16 cm. With these estimates, the calculated 
luminosity, vL v , is given as ~ 10 29 ergs s _1 which is well 
consistent with the observed values of 1.3 x 10 29 ergs s _1 
(0.1-2.4 keV) and 8.9 x 10 28 ergs s" 1 (0.5-10 keV). 

Although the general properties of the X-ray trail in 
PSR J2124— 3358 are not in contradiction with properties 
observed in other pulsars there are still ambiguities which 
are not completely resolved. First, one should notice that 
the trail is misaligned with the direction of the pulsar's 
proper motion. As reported by Gaensler, Jones, & Stap- 
pers (2002), the head of the H a bow shock is found to 
be highly asymmetric about the pulsar's velocity vector 
with the apparent nebular symmetry axis deviated from 
the velocity vector by ~ 30°. Even though the misalign- 
ment of the X-ray trail seems to agree with the asymmetry 
of the H a nebula, deeper observations by XMM-Newton 
and Chandra arc required in order to constrain the phys- 
ical properties of this interesting nebula in higher detail. 
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